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Abstract 

We study the Brownian momentum process, a model of heat con- 
duction, weakly coupled to heat baths. In two different settings of 
weak coupling to the heat baths, we study the non-equilibrium steady 
state and its proximity to the local equilibrium measure in terms of 
the strength of coupling. For three and four site systems, we obtain 
the two-point correlation function and show it is generically not mul- 
tilinear. 

Keywords: weak coupling limit, local equilibrium, Brownian momen- 
tum process, inclusion process, duality. 



1 Introduction 

In the study of non-equilibrium systems, exactly solvable models can serve 
as test-cases with which general statements about non-equilibrium, such as 
in f3], [11] can be tested. Recently, in [6], [7], P], we studied the Brownian 
momentum process (BMP) and showed that this models is exactly solvable 
via duality with a particle system, the symmetric inclusion process. In 
this paper, we look at the close-to-equilibrium states of the BMP. First, 
we consider a close-to-equilibrium scenario where the temperature of the 
right heat bath is close to the temperature of the left heat bath, and show 
that the distance between the local equilibrium measure and the true non- 
equilibrium steady state is of order at most the square of the temperature 
difference, in agreement with the theory of Mc Lennan ensembles, see |11] . 
Next, we consider a situation where the linear chain is coupled weakly to 
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heat baths to left and right ends (with fixed and different temperatures) , and 
study which equihbrium measure is selected in the limit where the coupling 
strength A tends to zero, as well as how far the true non-equilibrium steady 
state is from the local equilibrium measure for small coupling strengths. The 
temperature profile can be computed for all values of A and is only linear in 
the chain including the extra sites associated to the heat baths for A = 1, 
and linear if these sites are not included for all values of A > 0. Finally, 
we explicitly compute the two-point correlation for all A > for a three 
and four sites system and show that the multilinear ansatz of the two-point 
function introduced in [6], see also [3], [U fails for a system of four sites, 
except when A = 1. 

2 The model 

The Brownian momentum process on a linear chain {!,..., A'^} coupled at 
the left and right end to a heat bath is a Markov process {x{t) : t > 0} on the 
state space ^In = M'^^''^'-^^'. The configuration x{t) = Xi{t) : i G {1, . . . N} is 
interpreted as momenta associated to the sites i G {1, . . . , N}. The process is 
defined via its generator working on the core of smooth functions / : — )■ M 
which is given by 

TV 

L = \Bi + XBN + ^p{i,j)Lij (1) 

i,j 

with 

Lij — {xidj Xjdi) 

and where dj is shorthand for The underlying random walk transition 
rate p{i,j) is chosen to be symmetric and nearest neighbor, i.e., Pi^i+i = 
Pi+i,i = 1,^ £ {1; • • • 5 — 1}; p{h j) = Otherwise. Since Lij = Lj^i the 
symmetry of p{i,j) is no loss of generality. 

The boundary operators Bi,Bn model the contact with the heat baths, 
and are chosen to be Ornstein-Uhlenbeck generators corresponding to the 
temperatures of the left and right heat bath, i.e., 

Bi = ndl - xidi 

Bn = Tiid% — xnQn 

Finally, A > measures the strength of the coupling to the heat baths. 
The process with generator ([1]) is abbreviated as BMP\. 

If Tl = Tr = T, then, for all A > 0, the unique stationary measure of the 
process {x{t) : t > 0} is the product of Gaussian measures with mean zero 
and variance T. If / Tr there exists a unique stationary measure; the 
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so-called non-equilibrium steady state denoted by . The existence and 

uniqueness of the measure /u^^ j,^ follows from duality (see next section) . 
We will look at two different close-to-equilibrium scenarios: 

1. X = 1, Tr = Tl + e and e ^ 0, 

2. Tl / Tr, and A ^ 0. 

In both cases we look at the behavior of the measure /i^^ rp^, in case two, as 
A — )• 0, and in case one as e — )• 0. Since for A = 0, the system has infinitely 
many equilibrium measures, in the second case it is of interest to find out 
which of these measure is selected in the limit A — t- 0. Both in the first 
and second case, we want to understand how close the true non-equilibrium 
steady state is to the local equilibrium measure. 



3 Duality 

The BMP\ can be analyzed via duality. The dual process is an interacting 
particle system, the so-called symmetric inclusion process [8] , where particles 
are jumping on the lattice {0, 1, . . . , A'', A'"-!-!} and interacting by "inclusion" 
(i.e., particles at site i can attract particles at site j). The "extra sites" 
0, -|- 1 -associated to the heat baths- are absorbing. I.e., a dual particle 
configuration is a map 

e : {0,...,A^ + 1} 

specifying at each site the number of particles present at that site. The space 
of dual particle configurations is denoted by For ^ G 0^, ^^'-^ denotes 
the configuration obtained from ^ by removing a particle from i and putting 
it at J. 

The generator of the dual process then reads 

L,c/>(o = 2Aei[0(e''°)-c/>(O] + 

N-l 

+2Ae^[<^(e^'^+i)-c/>(0] (2) 

In words, this means particles at site i jump to j at rate 2p{i, j){2^j + 1). 
At the boundary site 1 (resp. N) particles can jump at rate 2A to the site 
(resp. N -\- 1) where they are absorbed. Absorbed particles do not interact 
with non-absorbed ones. The dual process is abbreviated as SIPx. The 
duality functions for duality between BMP\ and SIPx are independent of 
A and given by 

i?(e,x) = rl°T|-+^n 



=1 (2^^ 



1)!! 
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for ^ G Qff a dual particle configuration, and x ^ Qn- 
The duality relation then reads 

LD{^,x) = LaD{^,x) (3) 

where L works on x and on ^. By passing to the semigroup, from ([3]) we 
obtain the duality relation 

E,D{^,x{t))=ElD{m,x) (4) 

where E^^ is expectation in BMP\ starting from x G O^r, and E| is expec- 
tation in SIP\ starting from ^ G 

For G we denote |^| = Yl!i=o^ total number of particles in 

^. Since eventually all particles in a particle configuration ^ G 0^ will be 
absorbed, we have a unique stationary distribution fi^^ j,^ with 

f Ditx)fJ,^^,T^idx)= Yl TtTl,Fliat = ^) = k5o + l5N+i) (5) 

A:,«:fc+«=|,t| 

where ^(t = oo) denotes the final configuration when all particles are ab- 
sorbed and k5Q + l5N-\-i the configuration with k particles at and / particles 
at iV + 1. 



4 Temperature profile 

The local temperature at site i G {1, . . . , N} is defined as 

Ti = J xffi^^^T^idx) 

and by definition Tq = T^^Tn+i = Tr. We say that the temperature profile 
is linear in the lattice interval [K, L] if there exist a, 6 G M with Ti = ai + b, 
for all i G [K, L] . For the computation of the temperature profile we only 
need a single dual walker, which performs a continuous-time random walk 
with rates 2p(i,j) and absorption at rate 2A from the sites 1, A^. 
Indeed, using ([5]) we have 

Ti = TlFI (e(oo) = 6o) + (l - ¥l (e(oo) = 6o)) (6) 

From this expression, one obtains the following equations for the tem- 
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perature profile: 



N 



Y,P{hm = Ti-X{Tl-Ti) 

i=l 
N 



i=l 
N 

^pii,N)Ti = TN-XiTn-TN) (7) 

i=l 

The second equation expresses that the temperature profile is a harmonic 
function of the transition probabilities, whereas the first and third equation 
are boundary conditions. In the case A = 1 and p corresponding to the 
simple nearest neighbor random walk, the equation for Ti,i = 0,. . . ,N is 
the discrete Laplace equation, which gives a linear temperature profile in 
[0,A^ + 1]. 

Remark 4.1. In this paper we restrict to the symmetric nearest neighbor 
walk kernel p{i,j). The equations ([7]) hold for general symmetric p{i,j). 
However, in the cases where it is not translation-invariant and/or not near- 
est neighbor, the temperature profile will not be linear. 

We have the following theorem that follows immediately from the equa- 
tions dZD. 

Theorem 4.1. For all A > 0, the temperature profile is linear in [1, N] and 
is given by 

Ti = ai + b (8) 



1, . . ., with 



X{Tr-Tl) 



X{N -l) + 2 

Tl+Tr + X{NTl-Tr) 
X{N -l) + 2 

We can now look at different limiting cases: 

1. In the case A = 1 we recover the result from [6]: 

In this case (only) the temperature profile is linear in [0, + 1]. 

2. In the limit A — )• we obtain for all i G {1, . . . , N} 

hm r/^) = ^^±^ 
A^o * 2 
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3. In the limit A — t- oo we obtain Ti = Tl^Tj^ = T/j and the profile is 
linear in [1, A^], similar to a system with A = 1 and N — 2 sites. 



4. In the hmit — ^ oo, such that i/N — ^ r G [0, 1] fixed, 

lim Ti =Tl + r{TR - Tl) 

This means that the macroscopic profile is linear and does not depend 
on A. 

Remark 4.2. The expectation of the heat current in the steady state in the 
system is J = Tj+i — Tj. Heat conductivity k is defined via the equation J = 
kAT. From Theorem \4.1\ it follows that k = x(^j^\~^j^2 ^^^^/i is independent 
of the temperature (i.e. the system obeys the Fourier's law for all values of 

A>o;. 



5 The stationary measure for e — > 

We consider the first weak coupling setting, i.e, A = 1, Tr = + e. We will 
prove that up to corrections of order e^, the stationary measure is given by 
a product of Gaussian measures corresponding to the temperature profile, 
i.e., the local equilibrium measure. 

Let us denote this local equilibrium measure 

^Tl,Tr = ®f=iGT,{xi)dxi 

with Ti given by ([8]), 

GT(x) = -^exp(-xV2T) 

and iiti^,Tl+€ the true non-equilibrium steady state (with A = 1). Then we 
have the following result. 

Theorem 5.1. The true equilibrium measure and the local equilibrium mea- 
sure are at most order apart, i.e., there exists eo > such that for all 
^ G there exists a constant C = C{£,) < oo such that for all < e < eo 
we have 



j D{i^x)llTL,TL+e{dx)- j D{^,x)uTL,TL+eida 

Proof. For the local equilibrium measure we have 

N 



< aOe' (9) 



J D{C,xyT,,n+e{dx) = flTt (10) 
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expanding this up to order e we find, 



Start now from (jS]) and expand up to order e: 



' I 1 + ^ E l^i(^i^) = + ISn+i) I + 0{e') (11) 

fe,i:fc+«=|g| 

Upon identification of (jlOp and (jlip we see that we have to prove 
irf iCioo) = k6o + 16 N+i) = Eg(eoo(A^ + l)) 

k,l:k+l=\^\ 

The function (/>(^) := Eg(^oo(A^ + 1)) is the harmonic function for the dual 
process, i.e., 

L,d> = 



which satisfies the boundary conditions 

\ / ^ \ 

/ \i=l J 

Therefore, it suffices to show that 



^ 4=1 / \i=l / 



(N , ^, 



both satisfies 



and the boundary conditions (fT3|l . That satisfies the boundary conditions 
is immediately clear. The fact that tp is harmonic follows from explicit 



computation: 

= 26[^(e''°)-V'(0] 

7V-1 
i=l 



N + 1 

N-l 

+ Yl (2ei+i(2?i + 1)[-1] + 2^^(24^+1 + 1)[+1]) 

i=l 



i=l 



and since X]j=i ('^« ~ Ci+i) = S,i ~ we indeed have 

□ 



6 The case A ^ 

Next, we consider the second weak couphng setting, i.e., we fix ^ Tr 
and study the behavior of the measure /U^^ as a function of A. 

In this case, the local equilibrium measure is the product of Gaussian 
measures corresponding to the temperature profile dH]), i.e., we have to com- 
pare with where 

'^Tl,Tr = ®f=iGj.x{xi){dxi) 

where is given by ([8]). Denote 

m = I D{i,x) ^i^^^T^idx) (14) 

then (j) is the harmonic function of the dual generator satisfying the boundary 
conditions 

<^(e* = ^ + Mo + 15n+i) = HOMkSo + ISn+i) = 
On the other hand if we put 

m ■■= J D{^,^) 4,,Tnidx) = TtTi^HiT^^'^)^^ (15) 
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then we see immediately that ip satisfies the boundary conditions. 
We win now first prove 

Lemma 6.1. There exists Aq > such that for all ^ G Qj^ there exists 
A(^) > such that for all < X < Xq we have 

In particular, since there is only a finite number of dual particle configura- 
tions with total number of particles equal to K , we have, for all < X < Xq 

sup \{Ldij){0\<C{K)X^ 

for some C{K) > 
Proof. Compute 



Put Tr-Tn = Ti- Tl =: 7 



^ iV6.(2em + i)f%i-i 



Remember from Theorem 14.11 that Tj = Xai + 6, hence Tj — Tj+i = — Aa, 
with 

A(r^ - Tl) 



A(iV_i) + 2 

Tl+Tr + X{NTl-Tr) 
X{N-l) + 2 

We find 

Tr-Tl 
^ X{N-l) + 2 " 



and hence, from (fT6 
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We then see that the first order terms form a vanishing telescopic sum: 




and therefore; 



,W.(f)=4AVv,K)g(|±iii) 



□ 

Given this result, we will prove that the measures Vj^^ j,^ and /x^^ j,^ are at 
most order 0(Alog(l/A)) apart as A — )• 0. 

Theorem 6.1. Let (pjijj be the functions defined in (jl4|) and (|15p . then we 
have the following. There exists Aq > 0, such that for all ^ G il^ there is 
C(^) > 0, such that for all < \ < Xq 



i.e., in the limit A — t- 0, the equibrium measure corresponding to temperature 
{Ti + Tr)/2 is selected. 

Proof. We start with the following lemma 

Lemma 6.2. For all ^ £ $7^ a (dual) particle configuration, there exists 
c = c(^) > 0, a = a{(^) > such that for all A > 0, and for all t > 



\m-m\<c{o^^og\ 



(17) 



as a consequence. 




Proof. Using duality between BMPx and SIPx, and ([5]) 




¥.^D{i,xt) 4^,Tnidx) - / D{C,x) fi^^^Tnidx) 
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In order to see the last inequality, we remark that for a particle at posi- 
tions 1, A^, the probability to be absorbed at the next step is of order A, as 
the maximal rate to move to the other (non-absorbing) neighbor is at most 
2(1^1 + 1). □ 

Proof of Theorem l6.lt using Lemma 6.1, and duality between SIPx and 
BMPx, we have 



Combining with Lemma 6.2 we have 



Now optimize w.r.t. t by choosing t = (1/aA) log(a/A) □ 



7 The two point correlation functions in the hmit 

In this section we prove that for the two-point correlation function in the 
non-equilibrium steady state, the deviation from local equilibrium is of order 
A, which strengthens (fT7|l for ^ = 5i + 6j (i.e., we get rid of the log(l/A)- 
factor). In the appendix we give explicit expressions for the two-point func- 
tion of some finite systems, and show in particular that it is not multilinear 
for A / 1. 

Define for N} 

Yij = j {xfx'^j)^i^^^j,Jdx) 

and additionally Yoi = TlTj, Yi^^+i = TiTr. 

Denote by T the matrix with elements Tj^ = TiTj j and T^j = 'iTf 
\i i = j where Tj is the temperature profile of Theorem 14.11 

Theorem 7.1. There exists C > such that for all i,j G {1, • • • ,N} we 

have 

\Yij-Tij\<CX (19) 
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Proof. From the stationarity of fj,^^ we find that Y satisfies the fol- 
lowing system of linear equations for k^l G {1, • • • , -/V} 

= {-'^Yki + Yk_ii + Yk+ii + Yki-i + Yki+i) 

+ ^Ykk+i^ki + ^Yk-ikSiki - '^Yk-ikSik,i+i - -^Ykk+i^k,!-! 
+ X{TlTi - Yu)5ik + KnTk - Y^k)Sii 

+ X{TRTi-YNi)6Nk + X{TRTk-YNk)SNi (20) 
which has the form 

MY = D 

By explicit computation we obtain 

X :=M.T-D = 0(A2) (21) 
Prom this we will now derive that 



Y = T + 0(A). (22) 

Put 

||Y - T|| = ||M"^M(Y - T)|| = ||M"^X|| 
We will show that 

llM^^Xf < -^||Xf (23) 



which combined with (j2ip gives the desired result ([221 
To obtain (1231) consider 



< M-^X, M-^X > = < X, (M"^)'^M-^X > 
= <X,A-iX> 

with A := MM"^ Using the spectral decomposition of A, we get 
<X,A-iX> = ^^<X,ei><e„X> 

< — ^atII^ii' 

minj(A- ') 

where Aj are the eigenvalues of A with the corresponding eigenvectors 
So it suffices now to see that 

min(AS'^^) > cA^ 

i 

We have 

min(A,^^^) = inf <X,AX> 

i l|X||=l 



12 



The matrix M has the form M = K + AS and hence 



< X, AX >=< (K^ + AS^)X, (K^ + AS^)X > 



Therefore 



< X, AX > 
A2 



A^||S^X|p + 2A < S^X,K^X > +||K^X|p 
A2 



and so we obtain 



Um inf 



mmi(A. ') 
A2 



> 



Indeed, since M = K + AS is not singular, either K or S must not be sin- 
gular, therefore ||S^X||2 and ||K-^X||2 cannot be both zero. □ 

Remark 7.1. It follows from the correlation inequalities derived in that 
^ij > Tjj . Indeed, Tjj would be the correlation function if the dual walkers 
were walking independently, however, two dual walkers interact by inclusion 
(attraction), and this leads to a positive covariance. 
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9 Appendix 

Here we derive explicit expressions for the two point correlation function for 
systems with three and four sites. We start from the equations (|20p . 

Since Y^i is symmetric in k and I it suffices to consider k < I. the different 
cases are as follows; 



1. 1< k = 1 < N; {-2Ykk + 3Ykk-i + SYkt+i) = 

2. 1< k = 1-1< N-1; {-8Ykk+i+Yk^ik+i+Yk+ik+i+Ykk+Ykk+2) = 

3. l<k<l + l<N + l; {-4Yki + Yk-ii + Yk+n + Yki-i + Yki+i) = 

4. k = 1 = 1; (-2^11 + 3yio + 3^12) + \{TlTi - Y^) + \{TlTi - Y^) = 

5. k = 1 = N; {-2YNN + 'iYNN-l + 'iYNN+l) + KTRTN-YNN) + KTRTN- 

Ynn) = 

6. 1 = k = 1 - 1; (-8^12 + + 1^22 + Yii + yi3) + \{TlT2 - Y12) = 

7. k = 1-1 = N-1; {-SYn-in+Yn-2N+Ynn+Yn-in-i+Yn-in+i) + 
X{TrTn-i - Ynn-i) = 
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8. l = k<l + l<N + l; {-Wn + Yqi + Y21 + Yn-i + Yn+i) + KnT^ - 
Yii) = 

9. l=k<l + l= N + l; (-4yiiv + Foiv + Y2N + Y^-i + nAr+i) + 
\{TlTi - Y^m) + KTrTi - Y^n) = 

10. l<k<l + l = N + l; {-AYkN + Yk-iN + n+iJV + ifcJV-i + ifcJV+i) + 
\{TRn - YkN) = 



9.1 3 Sites System 



The equations for the two-point correlation function are of the form M.Y = D 
where 

/ -XTlTz - XTrTi \ 

-XTrT2 

-SXTlTi 
-XTlT2 




Y12 

Yl3 

Y22 
Y23 
\ Y33 ) 



and D = 



and the matrix M can be read from the previous equations as; 



1 -2(1 + A) 1 \ 

3 -(1 + A) 

11 -(7 + A) 1 
(1 + A) 3 

1 -(7 + A) 1 1 
VO 3 -2 3 0/ 

The expHcit solution is via inversion of M. The result for Y and the 
correlation functions Cjj = Yj^ — TiTj(l + 26ij) then reads as follows: 

_ 3 (r|(5 + 3A) + 2TlTr (5 + 9A + 2A^) + r| (5 + 23A + 24A^ + 4A=^)) 

4(1 + A)2(5 + A) 

rl (5 + 3A) + 2TlTr (5 + 4A + A^) + (5 + 13A + 2A2) 
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Yi3 = 



■ 22 



4(5 + 6A + A2) 

Tl (5 + 13A + 2A2) + r|(5 + 13A + 2X'^) + 2TlTr (5 + 9A + 12A2 + 2X^) 

4(l + A)2(5 + A) 

_ 3 {2TlTr (5 + 4A + A2) + r|(5 + 8A + A^) + (5 + 8A + A^)) 

4(5 + 6A + A2) 

r2(5 + 3A) + 2TlTr (5 + 4A + A^) + (5 + 13A + 2A2) 



23 



4(5 + 6A + A2) 
14 



3 (r2(5 + 3A) + 2TlTr (5 + 9A + 2\'^) + T|(5 + 23A + 2A\^ + 

I 33 — 



and 



4(1 + A)2(5 + A) 



?,{TL-TRy\ ^ {Tl-TrYX 



2(1 + A)2(5 + A)'^^^ 2(5 + 6A + A2) 



2(1 + A)2(5 + A)' 2(5 + 6A + A2) 

(r^-r^)2A 3(r^-TR)2A 

2(5 + 6A + A2)' 2(1 + A)2(5 + A) 

We see that for all A;, / 

Qu oc A(rL - r^)2 

and also Q^i > 0. 

One might be interested to see if the bi-linear ansatz introduced in [6] 
for the special case A = 1 is also valid here, i.e. 

Yij = a + bi + cj + dij 

Yii = A + Bi + Di^ (24) 

with the boundary conditions Yoj = T^Ti, Yj^^v+i = T^Tr. 

However, to check the validity of the ansatz we must calculate the cor- 
relation functions for a 4 sites system, since in 3 sites systems we have only 
6 correlation functions which are less than the 7 constants of the ansatz. 



9.2 4 Sites System 



Similar to the calculation for the 3 site system, we have M.Y = D where 



/ \ 

Yu 

Yl3 

Yu 
Y22 
Y23 

Y2A 
133 
134 

V 144 J 



and D 



/ \ 





- ATlTs 
-AT^ra 
-3ATlTi 
-XT2TR 

-XTrTs 
-SXT^Tr 
V -XTlT^ - XTiTr I 



and where the matrix M is given by 
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1 
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3 





-2 


3 








1 


-3- A 


1 





1 














1 


-7- A 


1 





1 

















-A 


3 



































1 





1 


-3- A 





1 























1 


1 


-7- A 


1 


























3 


-1-A 








1 


-2-2A 








1 












V 

The solution for Y is 



Y 
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6r|(12+A(14+3A))+6rLTR(24+A(76+5A(8+A)))+Tg(72+3A(172+3A(106+A(46+5A)))) 
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(6+A)(2+3A)(8+A(16+5A)) 
_ 2r|(l+A)(12+A(14+3A))+rj,ri;(48+A(152+A(128+A(44+5A))))+2T|(12+A(74+A(121+A(49+5A)))) 



Yi3 = 



14 



(6+A)(2+3A)(8+A(16+5A)) 
rg(l+A)(24+5A(4+A)(5+A))+T|.(24+A(76+A(41+4A)))+2rj.Tfl(24+A(76+A(109+A(47+5A)))) 

(6+A)(2+3A)(8+A(16+5A)) 
T|(24+5A(4+A)(5+A))+r2(24+5A(4+A)(5+A))+riTfl(48+A(152+A(314+3A(46+5A)))) 



^22 



(6+A)(2+3A)(8+A(16+5A)) 

_ 3(2Ti,Tfl(24+A(76+A(73+3A(9+A))))+T|(24+A(124+A(181+7A(10+A))))+T|(24+A(76+A(77+2A(12+A))))) 



23 



24 



(6+A)(2+3A)(8+A(16+5A)) 
2rg(12+A(5+2A)(10+A(8+A)))+2T|(12+A(5+2A)(10+A(8+A)))+rz,rfl(2+A)(24+A(64+A(50+7A))) 

(6+A)(2+3A)(8+A(16+5A)) 
_ T|(l+A)(24+5A(4+A)(5+A))+T|(24+A(76+A(41+ 4A)))+2Ti.Tfl(24+A(76+A(109+A(47+5A)))) 



Y33 

Y34 



(6+A)(2+3A)(8+A(16+5A)) 

3(2Ti.Tfl(24+A(76+A(73+3A(9+A))))+T|(24+A(124+A(181+7A(10+A))))+r|(24+A(76+A(77+2A(12+A))))) 

(6+A)(2+3A)(8+A(16+5A)) 
2r|(l+A)(12+A(14+3A))+TLTR(48+A(152+A(128+A(44+5A))))+2r|(12+A(74+A(121+A(49+5A)))) 

(6+A)(2+3A)(8+A(16+5A)) 
. _ 6T|(12+A(14+3A))+6rt,rR(24+A(76+5A(8+A)))+3r|(24+A(172+3A(106+A(46+5A)))) 
44 



(6+A)(2+3A)(8+A(16+5A)) 

and the corresponding correlation functions are 



p, _ 3(r^-Tfl)2A(24+A(50+13A)) ^ _ (T^-Tfl)2A(l+A)(24+A(50+13A)) 
~ (6+A)(2+3A)2{8+A{16+5A)) ' ^^12 - (6+A)(2+3A)2{8+A{16+5A)) 
_ 2(Tz.-Tfl)^A(l+A)(12+A(16+A)) _ 2(Ti,-rfi)^A(12+A(16+A)) 

'^IS ~ (6+A){2+3A)2(8+A(16+5A)) ' ^^14 — (6+A)(2+3A)2(8+A(16+5A)) 
f-^ _ 3(T^-Tfl)^A(2+A(4+A))(12+A(16+A)) ^ _ (Tt,-Tfl)2A(24+A(86+A(93+A(27+2A)))) 
^^22 - (6+A)(2+3A)^ (8+A(16+5A)) ' ^-^23 - (6+A){2+3A)2 (8+A(16+5A)) 

p, _ 2(Tj,-Tfl)2A(l+A)(12+A(16+A)) _ 3(T,,-rfl)2A(2+A(4+A))(12+A(16+A)) 

^24 - (6+A)(2+3A)2(8+A(16+5A)) ' ^^33 - (6+A)(2+3A)2(8+A(16+5A)) 
_ (Ti.-rR)^A(l+A)(24+A(50+13A)) _ 3(rj,-TR)2A(24+A(50+13A)) 

^34 - (6+A) (2+3A)^ (8+A(16+5A)) ' ^^4 - (6+A)(2+3A)^(8+A(16+5A)) 

We see once more that for all k,l 

Cki oc \{Tl - Tnf 



and Cjt/ > 0. 

Now we can directly check the validity of the bi-linear ansatz. Direct 
calculation shows that the diagonal part of the ansatz, i.e., Ya = A + Bi + 
Di^ is valid, but the non-diagonal part Yjj = a + bi + cj + dij is not. 
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If wc determine the coeficients a, b, c, d by fitting the bihnear ansatz to 
Yi2, Yi3, Y23, Y34, then we obtain 



(6 + A)(2 + 3A)(8 + A(16 + 5A)) ' 



which shows that the bihnear form can not hold for A ^ {0, 1}. Remark that 
also when A — t- 00 the deviation from the multilinear form vanishes, which 
is consistent with the intuition that this limit is the same as having A = 1 
in a smaller system obtained by removing the sites 1, N. 
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